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We investigate by means of Dynamical Mean-Field Theory the crossover from BCS supercon- 
ductivity to Bose-Einstein (BE) condensation of preformed pairs which occurs in the attractive 
Hubbard model by increasing the attraction strength. We follow the evolution of the two energy 
scales underlying the superconducting phenomenon, the gap Ao and the superfluid stiffness Ds, 
which controls the phase coherence. The BCS-BE crossover is clearly mirrored in a change in the 
hierarchy of these two scales, the smallest of the two controlling the critical temperature. In the 
whole intermediate-to-strong coupling region T c scales with Ds , while Tc is proportional to Ao only 
in the BCS regime. This evolution as a function of the interaction qualitatively resembles what 
happens in the cuprates when the doping is decreased towards the Mott insulator. 

This continuous change reflects also in the energetic balance at the superconducting transition. 
While, as it is well known, superconductivity is stabilized by a potential energy gain in the BCS 
regime, the strong-coupling superconductivity is made stable by a reduction of kinetic energy. Inter- 
estingly the intermediate-coupling region, where the maximum T c is achieved, behaves similarly to 
the strong-coupling regime, and its gain in kinetic energy is the largest as a function of the coupling. 
Since the integral of the optical conductivity is proportional to the kinetic energy, the above finding 
implies that the attractive Hubbard model can account qualitatively for the anomalous behavior of 
optical spectra around T c , where an increase of spectral weight is observed in under and optimally 
doped cuprates, while the overdoped samples have a more standard behavior. This qualitative 
agreement is lost in the normal phase, specifically at strong-coupling, calling for the inclusion of 
strong correlation effects in the theoretical description. 

PACS numbers: 71.10.Fd, 71.10.-w, 74.25.-q 
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I. INTRODUCTION 

High-temperature superconductors (HTSC) have at- 
tracted an unprecedented interest in the solid-state com- 
munity, yet the number of open issues largely exceeds the 
number of generally accepted points. Among the latter, 
we can count the crucial role of strong electron-electron 
repulsion. This immediately raises the question of the 
competition and coexistence between this repulsion and 
the attraction between quasiparticles responsible of the 
superconducting phenomenon. 

In this work, we make a step back, and we consider 
a simpler and more "phenomenological" question, that 
is the possibility to explain some features of the HTSC 
physics simply in term of a crossover between a conven- 
tional (BCS) superconductivity (characterizing the over- 
doped compounds) and a sort of strong-coupling Bose- 
Einstein (BE) superconductivity (in the underdoped 
region). 1,2 A number of properties of the cuprates can 
in fact be qualitatively described in this framework. No- 
table examples are the intermediate values (~ 10-^20 A) 
of the superconducting coherence length £o in optimally 
doped compounds 3 ^, and the pseudogap phenomenology 
in the underdoped compounds 5 , which can be ascribed to 
a "preformed pairs" phase, in which the superconducting 
order parameter has a finite amplitude, but it lacks phase 
coherence. The strongest connection with the physics of 



a BCS-BE crossover is perhaps the evolution with dop- 
ing of the relevant energy scales controlling the stability 
of the superconducting condensate, namely the gap Ao, 
which is proportional to the binding energy of the Cooper 
pairs, and the superfluid stiffness D$, which represents 
the energy cost for phase fluctuations. In the underdoped 
compound there is a clear experimental evidence^! that 
Ds < Ao and T c is proportional to the "weak link of 
the chain" Ds according to the so-called Uemura plot 8 , 
whereas for higher doping a more conventional direct pro- 
portionality of T c on Ao is recovered. 

However some obstacles arise trying to push fur- 
ther this line of thought. For instance, the uncon- 
ventional d— wave symmetry of the order parameter, 
with the related presence of quasiparticle excitations 
down to zero energy, is known to strongly affect the 
low-temperature thermodynamic properties. In partic- 
ular the nodal quasiparticles can dominate the low- 
temperature charg o 9 ! 10 and thermal transport even in the 
strong-coupling regime, invalidating a purely bosonic de- 
scription. 

It is also intuitively hard to reconcile the BCS-BE 
crossover scenario with the relevance of strong correlation 
effects, which are naturally larger and larger when the 
doping is reduced and the Mott insulator is approached. 
Then, in the BCS-BE scenario, the attraction should be 
stronger in the same region where the repulsion is maxi- 
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mum. This puzzling fact can be understood by interpret- 
ing the BCS-BE physics as relevant only for low-energy 
quasiparticles, while the high-energy physics has to be 
dominated by Coulomb repulsion. An explicit realization 
of a similar scenario has been obtained in models with 
orbital degeneracy and specific kind of interactional^*^, 
where the superconducting properties of the strongly 
renormalized quasiparticles share similarity with those 
of an effective attractive Hubbard modelii. 

Keeping in mind those limitations, we find it important 
to understand if, and to what extent, a simple attrac- 
tive picture for quasiparticles, modelized by an attractive 
Hubbard model, is able to reproduce some of the prop- 
erties of the cuprates. We will find that while some of 
the properties of the superconducting phase can be put 
in this framework, the normal state behavior necessar- 
ily requires the inclusion of strong repulsive correlations. 
The Hamiltonian of the attractive Hubbard model reads 



■h= i x '•:,'> 1 Y. 

<ij>a 



''■T 



where < ij > indicates that the first sum is restricted to 
nearest neighbors sites only, c\ a (c^) creates (destroys) 
an electron with spin a on the site i and = c\ a Ci„ is 
the number operator; t is the hopping amplitude and U 
is the Hubbard on-site attraction. As it can easily seen, 
this model reproduces in the extreme weak- and strong- 
coupling limit the BCS and BE regime respectively, and 
allows to move in the whole crossover region, simply by 
adjusting a single parameter, the ratio U/t. 

From the theoretical point of view, a stringent test of 
the relevance of the crossover in the attractive Hubbard 
model for the cuprates has been mainly limited by the 
lack of reliable non perturbative approaches able to fol- 
low the evolution from weak to strong coupling without 
a bias in some direction. In this work we use for this 
purpose the Dynamical Mean-Field Theory (DMFT)A£, 
a non perturbative approach which neglects the spatial 
correlations beyond the mean field level, but fully retains 
the local quantum dynamics, and becomes exact in the 
infinite coordination limit. The local nature of the at- 
traction and the non-perturbative nature of DMFT are 
expected to lead to reliable results and allow for a demo- 
cratic treatment of the different regimes and of the whole 
BCS-BE crossover. 

The attractive Hubbard model has been already inves- 
tigated with DMFT in the past, but the attention was 
focused mainly on the normal phase or on the determina- 
tion of the critical temperature for the superconducting 
transition, with the recent exception of Ref. Il4l where 
the BCS-BE crossover at T = is analyzed using iter- 
ated perturbation theory as an impurity solver, but the 
focus is mainly on spectral properties. More precisely, 
by explicitly avoiding superconducting solutions, a phase 
transition has been found both at fmit o 15 i 16 and at zero 
temperature^ between a metal and a "paired" phase, 
i.e., a collection of independent pairs without the super- 
conducting phase coherence. This paired phase represent 



the 'negative-?/' counterpart of the paramagnetic Mott 
insulator found for the repulsive Hubbard mode l 15 ' 17 . 

With the present work we complete the DMFT analysis 
of this model reported in the Refs. [illlllT7l with a care- 
ful investigation of the properties of the superconducting 
phase, certainly stable at low temperatures. Although 
the onset of the superconductivity smoothes the abrupt 
changes observed in the low-temperature metastable nor- 
mal phase, and the evolution of the superconducting 
phase as the coupling is increased is a smooth crossover, 
the way in which the energetics changes in this process is 
extremely interesting. We fully characterize the crossover 
from BCS to BE superconductivity, and we establish that 
the intermediate regime, where the maximum T c is ob- 
tained, shares the behavior of the strong-coupling, BE 
superconductivity. In particular, this means that the 
"optimal" superconductor is stabilized by a kinetic en- 
ergy gain, as it is expected in the bosonic regime, as 
opposed to the standard potential energy stabilization 
characteristic of weak-coupling superconductivity. 

The paper is organized as follows: Sec. II briefly de- 
scribes the DMFT for an s-wave superconductor, Sec. Ill 
is devoted to the evolution of the superfluid stiffness and 
the gap as a function of the interaction, Sec. IV presents 
the energetic balance at the superconducting transition, 
and the relation with optical measurements. Sec. V is 
dedicated to concluding remarks. 



II. METHOD 

In this section we will briefly introduce DMFT and 
some aspects related to the study of superconducting 
solutions. DMFT maps a quantum lattice model onto 
an impurity model whose hybridization function (usually 
called Weiss field in analogy with classical mean-field the- 
ories) is determined by means a self-consistent equation. 
The latter equation contains the information about the 
original lattice only through the non-interacting density 
of states. In our case Q is mapped onto an Anderson 
model with attractive coupling. In order to describe the 
superconducting phase, the bath presents superconduct- 
ing terms leading to an anomalous Weiss field 

Ham = ^ [ ei C L C ^ + Vi {c\ a d a + h.c.) 



+ Ai(cJ iC J T +h.c.) 



n, 



loc 



(2) 



where TLl 



Hoc = —U (not — |) ( n oi — |) — M n o is the on-site 
term, and the chemical potential fj, is adjusted to fix the 
particle density on the impurity site (in our calculations 
we always fix n = 0.75 as a g eneric density out of half- 
filling, as done in Refs. Ildl7h . 

From the impurity model we compute the normal and 
anomalous Green's functions, G(t) = — (Tcj(r)cj(O)) 
and F(t) = — (Tcj(r)cj(O)), which are used to build 
the matrix G{iu> n ) in Nambu-Gor'kov formalism. Anal- 
ogously one can define a matrix of Weiss fields whose 
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diagonal Q°(iu> n ) 1 and off-diagonal T a (iuj n ) 1 elements 
are related to the parameters appearing in @ by 

1=1 " I ' 

= +J2 V? J 2 - (3) 

i=l n "T e Z + 

The two above quantity also define the local self-energy 
matrix 

S(iw n ) = Q^iiuJn) - G _1 (iw„). (4) 

The self-consistency condition relates the Weiss 
field to the Green's function. We work with an 
infinite-coordination Bethe lattice with semicircular den- 
sity of states of half-bandwidth D, (i.e., N(e) = 
2/(ttD 2 ) V D 2 — e 2 ), for which the self-consistency reads 

Go (w n ) = iu n T° + [it 3 - t 2 T 3 G(iuj n )T3 (5) 

fj being the Pauli matrices. It is straightforward to check 
that Eqs. ©-© reduce automatically to their normal- 
state counterparts as soon A; = for each /. 

The inclusion of local quantum fluctuations rules out 
the possibility to analytically solve Eqs. 101-10, and re- 
quires numerical (or approximate) solutions of the Ander- 
son impurity model (J2J. Here we adopt Exact Diagonal- 
ization (ED) as the impurity solver. Thus, we discretize 
the Anderson model, by truncating the sums in Eqs. (0) 
and 10 to a small finite number of levels N s . It has 
been shown that extremely small values of N s provide 
really good results for thermodynamic properties and re- 
liable results for spectral functions^ Here we use both 
the Lanczos algorithm (at zero temperature) and the fi- 
nite temperature algorithm in its simplest version, which 
requires the full spectrum of the Hamiltonian matrix. To 
obtain the full spectrum of the Hamiltonian, needed to 
compute the finite temperature properties, we are forced 
to a rather small value of n s , namely 5, whereas in the 
zero temperature case the Lanczos algorithm allow us to 
deal with larger systems, up to n s ~ 10 also in the super- 
conducting phase, where the number of particles is not 
conserved, making the size of the Hilbert space larger. 



III. HIERARCHY OF ENERGY SCALES 

As we anticipated in the introduction, a landmark of 
the HTSC phase-diagram is the crossing of the energy 
scales relevant for superconductivity: the superconduct- 
ing gap Ao and the superfluid stiffness D$, which repre- 
sent respectively the energetic cost to form a Cooper pair 
and to break the phase coherence of the superconducting 
phase. 

DMFT allows for a straightforward calculation of the 
order parameter Aq at zero temperature, defined through 



the local anomalous Green's function 

A = -£/<Tc T (0) Ci (0)) = -U/f3j2n^n)- (6) 

n 

We note that the values for Ao obtained with this stan- 
dard "mean-field" definition coincides almost exactly 
with the anomalous part of the self-energy £i2(iw„) at 
large u n . In our approach £12 is frequency dependent, 
with a small reduction at small Matsubara frequencies 
with respect to the large uj n value. In most cases we find 
that the above reduction is quite small and £12 — Ao. 

On the other hand, the superfluid stiffness D$, is de- 
fined in terms of the static limit of the electromagnetic 
response function as 

D s = D dia - Xjj (ci^0,fl = 0). (7) 

The diamagnetic term D^ia is given by Ddia = 
-{E Hn ) = -2/pJ2 u JdeeN(e)G{e,LJ n )&). The para- 
magnetic term, which measures the normal component 
Dm is defined as the sum over all the directions^ of the 
transverse part of the paramagnetic kernel in the static 
limit (i.e., X jM - 0,Q = 0) = V„ C' ;<1 -> 

O.M = 0) being x^(r,r) = (T T j Q (r, r)j a (0, 0))). 

In principle the evaluation of XjjiQ — * 0, O = 0) would 
require the calculation of the corresponding four-field cor- 
relation function, but a remarkable simplification occurs 
in the infinite coordination limit where the DMFT is ex- 
act, since all the vertex corrections to the e.m. kernel 
vanish 1 -. The evaluation of D$ requires therefore only 
the dressing of the non-interacting Green's function in 
the electromagnetic kernel with the self-energy. As a re- 
sults, the DMFT expression for \jj readsi&SI 

" n J 

x [G(e, cj„)G* (e, w n ) + F(e, u> n )F{e, w n )] 

where V(e) = (At 2 — e 2 )/3 is the square current ver- 
tex for the Bethe lattic o 21 ' 22 i 2 ? and G(e,u n ), F(e,w n ) 
are the normal and anomalous lattice Green function 
respectively 2 ^. 

Thanks to the specific form of the current vertex, it 
is possible to write a more compact expression for D$- 
More specifically by exploiting the relation —eN(e) = 
d e [V(e)N(e)] and then transferring the energy-derivative 
on the Green functions one finally gets 

D S = I £ / de N ( e ) V ( e ) F ( e > Wn) (8) 

The above expression only contains the anomalous Green 
functions, and makes it apparent the vanishing of Ds at 
T = T c . 

Before comparing the evolution of Ao and Ds, we 
briefly discuss the temperature behavior of the super- 
fluid and normal densities. In particular, the evaluation 
of D/v and Ds reveals that the T = value of Dm is basi- 
cally negligible, as predicted by BCS mean-field, not only 
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FIG. 1: (Color online) Temperature dependence of the su- 
perfluid stiffness Ds and of the normal component Dn at 
U — 2.0D (just before the maximum T c ) compared with the 
diamagnetic term — (Ekin)- 



FIG. 2: (Color online) Superconducting energy scales (in 
unit of D) (Ao and Ds have been normalized to be directly 
comparable with T c in the asymptotic regimes). The two 
energies cross at intermediate coupling at U ~D, well before 
the maximum T c dome (T c is taken from Ref. uff\ 



in the weak-coupling regime, but also for sizable interac- 
tion, as shown in Fig. 2] for U = 2D, where D^/Ds at 
T = turns out to be less than 0.01. This means that it 
is possible to identify D$(T = 0) with — (Ekin) in a wide 
range of couplings. It is clear that this results depends on 
the neglect of small-momentum collective modes, which 
can deplete the condensate 25 . It is worth noting that the 
presence of a frequency-dependent self-energy and of an 
incoherent part of the Green's function does not lead to 
a reduction to a depletion of Ds, differently from what 
happens, for instance, for impurity scattering. 

We now come to the evolution of the energy scales in- 
herent to the superconducting phase as a function of the 
interaction. Our results are summarized by Fig. [21 where 
the superconducting gap Ao, the superfluid stiffness Ds 
and the critical temperature T c are plotted as a function 
of the ratio U/D. Ao is smallest in the weak-coupling 
regime, and the superfluid stiffness becomes the lowest 
scale in strong-coupling. It is natural that the system 
becomes superconductor only when the pairs are formed 
and they have phase coherence, namely when the tem- 
perature is lower than both the gap scale and the super- 
fluid stiffness scale. Therefore the critical temperature is 
substantially determined by the smaller scale at each in- 
teraction. Thus the critical temperature is proportional 
to Ao in the weak-coupling regime, as predicted by BCS 
mean-field, while in the strong-coupling limit we recover 
T c ~ D 2 /W ~ (E km )/2 ~ As/2, as predicted by the 
mean-field solution in the limit of an hard core boson 
system with Heisenberg coupling J = D 2 /(dU)r^ 

The crossing between the two energy scales, which is 
characteristic of the BCS-BE crossover occurs at interme- 
diate coupling for U slightly smaller than maximum T c 
value. Therefore only in the proximity of the weak cou- 
pling regime the smallest energy scale, which ultimately 
determines T c , is the binding energy of the Cooper pairs, 
whereas already at the optimal U (and of course in all the 
strong coupling regime) the weak link for the stability of 



the superconducting phase is the superfluid stiffness Ds 
(hence, the onset of the superconductivity is controlled 
by the phase-coherence). In other words, the optimal su- 
perconductivity is achieved in a regime where the physics 
is already that of strong-coupling, with pair formation 
occurring at a temperature higher than T c . We notice 
that the optimal interaction value is noticeably larger 
than the attraction needed to form a bound state in the 
low-density limit Ub ~ fli». 

Further insight on the influence of the crossing of 
the energy scales on T c can be highlighted by plotting 
the critical temperature versus Ds- In this way a sort 
of "Uemura-plot" relation, like that characterizing the 
physics of the underdoped cuprates^, is found in the at- 
tractive Hubbard model for U > 2 AD, where T c is a 
decreasing function of U (see Fig. |3J). This result sup- 
ports the identification of the superconductivity with a 
condensation of preformed bosons already at a moderate 
values of the interaction. 

On the other hand, the attractive Hubbard model, at 
least in our DMFT framework, cannot reproduce the de- 
viations from the Uemura plot recently observed in the 
cuprates2Zi2&. More specifically an evident "re-entrance" 
of the curve T c (Dg) has been observed in many over- 
doped compounds, and its shape appears to be strongly 
material dependent. In our calculation at small U, where 
Ds — — (Ekin), the superfluid stiffness is a monotonically 
decreasing function of the interaction, and does not show 
the experimentally observed re-entrance. 

In order to push a little bit further the comparison 
with the most recent experimental data, it is useful to 
recast our DMFT results for the "Uemura plot" rescal- 
ing T c with the superconducting gap $o in the single 
particle spectrum, as done by Tallon et aS^, in order to 
single out as much as possible the nonstandard behav- 
ior of T c . The experimental T c /$o is found sub-linear in 
Ds in the underdoped regime, while it is slightly super- 
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FIG. 3: (Color online) Uemura plot in the attractive Hubbard 
model. The direct proportionality between Ds and T c in the 
strong-coupling limit is evident on the left side of the figure, 
corresponding to U > 2.4D. 



linear at higher doping (with a not universal tail in most 
overdoped compounds). We notice that, beyond the BCS 
regime $o can be different from the order parameter A 
that we have introduced before. Here $o is evaluated di- 
rectly as the gap in the density of states. The difference 
between the two quantities varies up to 30%. 

Our results are summarized in Fig. 0] where one can 
note the presence of a wide region of a sub-linear behav- 
ior of T c /4>o on the left-side of the figure, which corre- 
sponds to the intermediate-to-strong coupling regime in 
fair agreement with Ref. IH- Also a partial of a super- 
linear behavior appears on the right-hand side of our plot. 
The observed behavior of T c j $ is therefore captured by 
our simple BCS-BE crossover, without invoking more in- 
volved explanations^. 

The evolution of energy scales we have discussed in 
this section is reflected in an increased role of phase fluc- 
tuations as the coupling grows. As we have mentioned 
in the introduction, those fluctuations may have in prin- 
ciple relevant effects on the superconducting phase, in 
which finite dimensionality effects beyond DMFT can be 
important, possibly destroying the superconducting or- 
dering. In DMFT, indeed, phase fluctuations at large q 
are properly taken into account, while the small q col- 
lective modes are neglected. In Ref. it is argued that 
superconductivity is destroyed by phase fluctuations for 
U larger than a critical value due to the small charge 
compressibility k — dn/dfi oc l/U derived within a 
phase-only effective theory based on the atomic limit, 
. Since k measures the inertia of the system against the 
dynamic phase fluctuations, a small value would permit 
large zero-point fluctuations, eventually destroying phase 
coherence. On the other hand, an alternative derivation 
of a phase-only action within one-loop expansion gives 
instead k ~ 2U/D 2 , therefore increasing with U—. 

Even if DMFT neglects finite-dimensionality effects 
and cannot therefore treat the physics of the Goldstone 
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FIG. 4: (Color online) Rescaled "Uemura-plot" (following 
Ref. l2Si) for the attractive Hubbard model. We have plotted 
the T c divided by the gap $o as a function of the superfluid 
density. It is evident the sub-linear behavior of T c /<3?o in the 
intermediate-to-strongcoupling regime, which mirror rather 
well the data of Ref. [2|| for several underdoped or optimally 
doped cuprates. 
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FIG. 5: (Color online) Charge compressibility in the super- 
conducting phase at T « T c . The DMFT data are com- 
pared respectively with the strong-coupling behavior from 
the phase-only action of Ref. [2^ and with the noninteract- 
ing value. 



modes associated with the phase-fluctuations, we can use 
this approach to evaluate the coefficients appearing in 
the phase-only effective theories. The lack of bias to- 
wards weak or strong coupling of DMFT should allow us 
to discriminate between the two derivations. The results 
shown in Fig. [51 clearly display that k grows line arly 
with U at strong coupling in agreement with Ref. |23 
the cost of dynamic phase fluctuation at q ~ tends to 
increase in the large U limit, leaving the sup erconducting 
phase stable, in contradiction with Ref. l29l A more com- 
plete description of the phase-fluctuations would require 
the calculation of all the phase-only theory coefficient, 
including the anharmonic ones. This issue is beyond the 
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aim of the present work. 



IV. ENERGETIC BALANCE OF THE 
SUPERCONDUCTING PHASE AND OPTICAL 
SUM RULE 

The characterization of the superconducting phase can 
be made more concrete by studying the temperature be- 
havior of the kinetic (Eki n ) and the potential (E pot ) en- 
ergies. From one side this analysis allows to evaluate sep- 
arately all the energetic contributions which characterize 
the mechanism for the stabilization of the superconduct- 
ing long range order across the BCS-BE crossover. On 
the other hand, since in a lattice system the frequency 
integral of the optical conductivity cr(u>) can be related to 
the kinetic energy of the carriers^, our calculations allow 
for a direct comparison with the optical measurements on 
the cuprates. 
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A. Kinetic and Potential Energy behavior 

The calculation of both potential and kinetic energies 
is quite straightforward in DMFT, since it only involves 
local quantities. This is evident for the potential energy 
because by definition it is proportional to the local den- 
sity of double occupancies n c i = (J^ ^it n «'l) 

E pot = ~Un d . (9) 

As far as E^m is concerned, exploiting the simplified 
form of the self consistency equation for the Bethe lattice 
we can derive the following expression 

Ekin = t 2 ry^[G(w n ) G(iu) n ) + G*(iu n ) G*(iu n ) 

iuj n 

- 2F(iu n ) F(iu n )], (10) 

which shows that also the kinetic energy can be expressed 
only in terms of the local (normal and anomalous) Green 
function. 

In Fig. we report the kinetic energy and potential 
energies for U = 0.8D, 2 AD and 3.6D, chosen as rep- 
resentative of the BCS, the intermediate and the BE 
regimes, as a function of the temperature. The on- 
set of the superconductivity is always marked by an 
abrupt change in the concavity of Eki n {T) e E pot (T). 
The energetic balance of superconductivity displays in- 
stead clear differences between the various regimes. At 
U = Q.8D the onset of superconductivity is accompa- 
nied by a slight loss in the kinetic energy (A£?^n — 
Ekin(0) ~ E kin (T c ) = +0. 003D), and potential energy 
gain (AEp t = E pot (T = 0) - E pot {T c ) = -0.008D), as 
in standard BCS theory. On the other hand, both at 
U = 2AD and U = 6AD the superconducting phase is 
characterized by a lower value of E^n an d a loss of po- 
tential energy is observed below T c . 



FIG. 6: (Color online) Low-temperature behavior of the ki- 
netic (upper panels) and the potential (lower panels) energies. 
The critical temperature is marked by arrows. Notice that the 
kinetic energy is normalized by the half-bandwidth, and the 
potential energy by the attraction strength U 

Such changes in the energetic balance at the supercon- 
ducting transition clearly highlight the different mecha- 
nisms stabilizing superconductivity in the two regimes. 
In the BCS limit superconductivity coincides with pair 
formation, which determines a gain in potential energy 
and a consequent loss in kinetic energy. In the opposite 
BE regime, the electrons are paired at a high tempera- 
ture of order U, but a true long-range superconducting 
order can take place only when phase coherence estab- 
lishes between pairs. Therefore T c is associated to a gain 
in kinetic energy, while a small fraction of the potential 
energy gained with pair formation is lost at T c . 

More interestingly, our calculations indicate that al- 
ready at intermediate U the "strong-coupling" mecha- 
nism is taking place, and superconductivity is associ- 
ated with a consistent kinetic-energy gain, such as the 
one observed in the optimally doped cuprates (see next 
sub-section). The kinetic energy gain turns out to be 
maximum at intermediate U/D, where the highest T c is 
achieved. The overall picture is drawn in Fig. [7] where 
we have reported the variation of the kinetic AEkin , 
the potential AE pot and the total energy AE to t between 
T = T c and T = for a number of interaction values. 
From these data one can easily see that (i) the typical 
BCS feature of an increasing E^in below T c is rapidly 
lost (approximatively for U ~ D ~ Ub), in agreement 
with the discussion of the previous section; (ii) an inter- 
mediate region exists, in which the onset of the SC phase 
is associated with a gain of both potential and kinetic 
energy; (in) the BE regime is effectively reached for U 
larger than X.8D (hence before reaching the maximal crit- 
ical temperature), where superconductivity is stabilized 
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by kinetic energy. 

It is also worth noting in Fig. [TJthat the variation of the 
total energy AEtot has approximatively the same dome- 
shape behavior of the T c dependence on U. It is tempting 
then, to consider this value as a measure of the conden- 
sation energy E c = E^ ot (T = 0) - Eg t (T = 0) of the su- 
perconducting phase. This identification works well only 
if one can assume a small variation of E^ ot (T) between 
T = T c and T = 0, an assumption which holds both at 
small (due to the small value of T c ) and at intermediate- 
large U (where both Ekin and E po t are almost constant 
in temperature above T c ). However this is not the case 
for U ~ 1t2D, where a strong temperature depen- 
dence of Ekin(T) (and of E pot (T)) in the normal phase 
is found, because of the presence of a strongly renor- 
malized quasiparticle excitation at the Fermi leveliS^i. 
Within DMFT the situation is even more involved for 
lid ~ 2.2D < U < 2.9D ~ U c % where a coexistence 
between metallic and insulating solution is found and 
the temperature dependence of E k i n and E pot in the 
normal phase below T c is subject to extremely abrupt 
changesiS^i. However this is a peculiar feature of the 
DMFT treatment of the Hubbard model which could not 
be so relevant when comparing our result with the exper- 
imental findings, also because these features of the Hub- 
bard physics occur at temperatures much smaller than 



B. Comparison with the optical measurements 



The integral over all the frequencies of the real part 
of the optical conductivity in a lattice model is re- 
lated to the second order derivative of the free particle 




FIG. 7: (Color online) Kinetic, potential and total energy 
variation in the superconducting region for different values of 
the pairing interaction U 
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and the effect of the interactions is hidden in the val- 
ues of the momentum distribution function n^ a . Here 
(Ekin) = jf Ek, ff e k«-k : <r , a the lattice spacing, d the 
dimension of the system considered, and V the unit-cell 
volume. The second equality, which holds for nearest- 
neighbor hopping only, implies that the kinetic energy 
data in Figs. El and [7| can be compared directly with the 
the optical spectral weight of the conduction band in the 
cuprates. 

The experimental results for the low-energy behavior 
of the spectral weight behavior in the HTS have been 
the object of a recent intense debate, and consensus has 
been reached about a few general points. In all the 
cuprates a sensible enhancement of low-energy spectral 
weight W(T) = L c diua(uj,T) (Tic being a frequency 
cut-off of the order of 1 eV, which is the plasma fre- 
quency for those materials), is observed when lowering 
the temperatur o 30 i 33 ' 34 i 35 . In the superconducting phase, 
underdoped (UD) and the optimally (Opt) doped sam- 
ples behave differently from the overdoped (OD) ones. 
More precisely, in the UD-Opt compounds below T c 
an upward bump is observed with respect to the nor- 
mal phase behavior, in contrast with the decreasing of 
W(Q, C ,T) observed in OD samples^. 

The gain of kinetic energy per Cu atom AE k i n — 
E kin (T c ) - E km (0) in UD-Opt BSCCO can be estimated 
using (jTTJ) as AE kin = -(0.5 -h 1) me V 30 i 33 i 36 . Quite 
noticeably, this value of AE k i n is much higher than 
the condensation energy for the same material 3 ^, esti- 
mated as E c ~ 0.1 meV on the basis of the specific heat 
measurements 3 ^. 

According to our DMFT results, the bump of 
W(£lc,T) observed below T c , can be qualitatively un- 
derstood in the framework of the BCS-BE crossover. 
More specifically the hypothesis of an intermediate-to 
strong coupling description of the superconductivity in 
the UD-Opt region appears to fit well with the observed 
enhancement of the low-frequency spectral weight below 
T c , whereas the disappearance of the upward bump in 
the OD cuprates is perfectly compatible with a weak- 
coupling superconductivity. It is interesting therefore to 
check to what extent is possible to push forward such an 
analogy with the properties of the real systems. 

We can try to establish a comparison with the experi- 
mental data by choosing the value of the semi-bandwidth 
D of the attractive Hubbard model in order to reproduce 
the maximum value T™ ax (a T™ ax of 90 - 100 K is ob- 
tained with D - lOOOif - lOOmeU). The theoretical 
values for |AE^j n | and |Ai? tot | (that we take as an esti- 
mate of E c ) are respectively 4-^8 meV and 2-^4 meV, 
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larger by less than one order of magnitude than the ex- 
perimental values in the underdoped HTSC. It is anyway 
remarkable that the DMFT of the attractive Hubbard 
model correctly predicts the qualitative trend of a con- 
densation energy which is only a fraction of the kinetic 
energy gain. 

As soon as the temperature exceeds T c and we enter 
the normal phase, our calculation ceases to properly de- 
scribe the experimental results. The T 2 behavior of E^ in 
is in fact found only in the BCS region, as shown by the 
U = 0.8D data of Fig. H3 and it completely disappears 
in the intermediate-to-strong coupling regime: no rele- 
vant variation of Eki n (T) is found at U — 2.4, 6. AD in 
a wide temperature range above T c . The inadequacy of 
the present approach to describe the experimental find- 
ings in the normal phase has different origins. A first 
reason is the freezing of spatial fluctuations characteris- 
tic of single-site DMFT. It is indeed reasonable to expect 
a role of short-range fluctuations at least at strong cou- 
pling. In this limit the attractive Hubbard model can be 
mapped in an effective "pseudospin" hamiltonianl, and 
the kinetic energy is given basically by nearest-neighbor 
correlations between the pseudospin operators, i.e., the 
local pairs and the empty sites. The temperature depen- 
dence of these correlations is poorly described by DMFT 
at large U/D. The relevance of this effect can tested 
using cluster extensions of DMFT^&Si, where the short- 
range spatial correlations are included. 

There is however a more basic reason for the in- 
adequacy of the attractive Hubbard model to capture 
the temperature behavior of the normal phase of the 
cuprates, which is the inability to properly describe the 
approach to the Mott insulator as the doping is reduced. 
The underdoped region is indeed almost universally be- 
lieved to be dominated by the physics of a doped Mott 
insulator, with strongly renormalized quasiparticles. On 
the other hand, the strong-coupling region of our BCS- 
BE framework is certainly a "correlated" regime, with 
renormalized quasiparticles, but the interaction is sim- 
ply the attractive one, which does not lead to the Mott 
insulating state. Indeed in a recent work^i it has been 
explicitly shown that the physics of a doped Mott insu- 
lating system, with its associated band-narrowing, may 
represent the natural explanation for the strong temper- 
ature dependence of the spectral observed in the normal 
phase of the cuprates. 



V. CONCLUSIONS 

In this paper we have performed a detailed investiga- 
tion of the physics of superconducting phase in the at- 
tractive Hubbard model both at zero and at finite tem- 
perature. The approach used here, namely the Dynam- 
ical Mean Field theory, is completely non-perturbative 
allowing for a treatment of the superconducting phase 
properties, which is not tied in principle either to the 
weak-coupling (BCS) or to the strong-coupling (BE) 



regime. 

In particular we have investigated the behavior of the 
energy scales relevant for the superconductivity as a func- 
tion of the pairing interaction U. We have found that the 
evolution of the superfluid density Dg and the supercon- 
ducting gap Ao displays a clear crossing in the intermedi- 
ate coupling regime, which is reminiscent of what is found 
experimentally in the phase-diagram of HTSC for differ- 
ent doping levels. In the weak-coupling region the BCS 
picture works well and the superconductivity is controlled 
by the binding energy of the Cooper pairs Ao; on the con- 
trary when the pairing interaction is high the supercon- 
ductivity becomes essentially a phenomenon of superflu- 
idity of preformed local pairs. This is clearly witnessed by 
the direct proportionality found in this regime between 
T c and D$, and showed in a sort of "Uemura-plot" for 
the attractive Hubbard model. We also consider explic- 
itly the problem of superconductivity suppression due to 
phase fluctuations in the hydrodynamic regime, which 
is usually neglected in a DMFT framework. Our DMFT 
calculation of the compressibility indicates that, contrary 
to the claim of RefiSS, the phase-fluctuation effects do not 
destroy superconductivity, even for extremely large val- 
ues of U. 

The analysis of the superconducting phase properties 
in the attractive Hubbard model has been then further 
enriched by considering the energetic stabilization of the 
superconductivity and the f sum-rule of the model. Our 
results clearly show a remarkable difference in the en- 
ergetic balance responsible for the stabilization of su- 
perconducting order when moving along the BCS-BE 
crossover. While in the weak-coupling regime the onset 
of superconductivity is associated to a gain of potential 
energy and a slight loss in kinetic energy as it is expected 
in a BCS picture, already in the intermediate coupling at 
U ~ D and specifically in proximity of the maximum of 
the T c dome, the situation is completely reversed. Here, 
and for higher values of U, the stabilization of the super- 
conductivity is due to a marked reduction of Ekm , while a 
smaller (but sizable) loss of potential energy is observed. 
Therefore the onset of superconductivity is clearly dis- 
tinct from Cooper pair formation already for moderate 
values of the pairing interaction. 

The energetic balance at the superconducting transi- 
tion finds a direct experimental counterpart in the optical 
measurements, since the integral of the optical conductiv- 
ity is proportional to the kinetic energy itself. Our results 
can be rephrased in this terms, so that in weak-coupling 
a slight reduction of the optical sumrule is found in our 
calculation at the superconducting transition, while a rel- 
evant enhancement of the sumrule is observed for T <T C 
in the whole intcrmediate-to-strong coupling region. This 
scenario is actually realized in optical measurements in 
the HTSC cuprates as a function of doping, with over- 
doped samples behaving more or less as standard BCS 
superconductors, and underdoped one, characterized by 
a gain in kinetic energja^aSl. 

A more quantitative comparison with the experiments 
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reveals that the variations of both kinetic energy and con- 
densation energy between T = and T c are reasonably 
well reproduced by our model, which also correctly pre- 
dicts that a sizable fraction of the kinetic energy gain is 
canceled by the potential energy loss, leading to a con- 
densation energy significantly smaller than the kinetic 
energy gain, as it indeed happens^ in the HTSC. It 
must be noticed that our calculation tends to overesti- 
mate both kinetic energy gain and condensation energy. 
However, the main limitation of the attractive Hubbard 
model in this context is the inability to capture the T 2 
enhancement of the spectral weight in the normal phase 
above T c . The main reason for this disagreement is likely 
the complete neglect of the strong repulsive interactions 
characterizing the underdoped cuprates close to the Mott 
state. 
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FIG. 8: (Color online) Same quantities of Fig. 2 plotted as 
a function of D/U. Both Ds and Ao are normalized by U. 
Now the crossing between the superconducting gap Ao and 
the superfluid stiffness at D ~ U, occurs in the proximity of 
the maximum value of T c . 

An important outcome of our analysis is that the op- 
timal critical temperature is obtained for an interaction 
strength well inside the BE region. This finding is in 
good agreement with previous calculations of the pseu- 
dogap temperature T* extracted from the spin suscep- 
tibility and specific heat, where the T* is significantly 
larger than T c at optimal interaction^, but it is in con- 
trast with experiments in Bismuth and Yttrium based 
cuprates, where T* tends to vanish close to optimal dop- 
ing. It is interesting to notice that this discrepancy seems 
to disappear if we fix the attraction strength U and fol- 
low the BCS-BE crossover by varying the half-bandwidth 
D, and measuring energies in units of U. Such a rescal- 
ing is meant to roughly describe a situation in which 



the crossover from BCS to BE is due to a shrinking of 
the coherent band due to strong repulsive correlation ef- 
fects which are stronger and stronger as the Mott insula- 
tor is approached, while the attraction, which we might 
think to arise from antiferromagnetic superexchange for 
the sake of definiteness, is basically unrenormalized in 
the same process. 

As shown in Fig. [SJ the maximum critical tempera- 
ture now occurs for U ~ D, at the boundary of the BCS 
region, where the pseudogap is small. However this dif- 
ferent perspective, in which the optimal T c moves closer 
to the BCS region is in contradiction with the experimen- 
tal evidence of a kinetic-energy driven superconductivity 
around optimal doping, which is qualitatively obtained 
in the BE regime. 

As a summary of our results, we can draw some fi- 
nal considerations on the the connection of the BCS-BE 
crossover scenario and the physics of the HTSC: several 
qualitative feature of the HTSC properties (as the cross- 
ing of the hierarchy of Ds and Ao, the Uemura plot ,. . . ) 
can be actually captured within a purely attractive de- 
scription. However, as a general trend, it provides too 
high estimation of the values of many thermodynamic 
quantities as the kinetic energy variation below T c and 
the condensation energy. At the same time the purely 
attractive description fails in reproducing the magnitude 
of the temperature dependence of the optical spectral 
weight above T c , calling for the inclusion of strong corre- 
lation effect and suggesting that a closer agreements with 
the physics of the cuprates may be obtained by applying 
the attractive Hubbard description, or more generally the 
BCS-BE crossover picture, only to the narrow quasipar- 
ticle excitations characteristic of a doped Mott insulator. 
Models in which explicit attractive and repulsive inter- 
action are presenljAiii^ are natural candidates to explore 
this scenario. 

After completion of this work, we became aware of 
recent Cellular DMFT 39 results by Kyung, Georges and 
Tremblay (cond-mat/0508645), where it is proposed that 
the inclusion of short range correlations beyond single- 
site DMFT makes the evolution on the normal state con- 
tinuous, and the zero-temperature energetic balance of 
the superconducting transition is studied. 
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